Wiecej, niz Gra w Zycie 
(SmoothLife) 


Modelowanie komputerowe 


Maciej Matyka 


Plan 


ο Automaty komórkowe ogólnie 

ο Sąsiedztwo 

ο Plamy 

e Gra w życie 

Piasek 

e Gaz HPP 

ο Więcej, niż gra w życie (Larger than life) 

ο SmoothLife 

ο Elementarne Automaty Komórkowe (Wolfram) 


Automat komörkowy 


Proste zasady lokalne -> globalne efekty w uktadzie 


Czym jest automat komórkowy? 


е Model dyskretny w czasie 

ο Model dyskretny w przestrzeni (1D, 2D itd.) 
ο Stan układu = wartości w komórkach 

ο Synchroniczna zmiana stanu komórek 

ο Zmiana stanu komórki zgodnie z regułą 


e Reguła - funkcja stanu komórki i otoczenia 


Historia 


ο Stanistaw Ulam 
https://en.wikipedia.org/wiki/Stanislaw Ulam 


е John von Neumann - 1948 - 
https://en.wikipedia.org/wiki/John von Neumann 


ο John Conway - 1970 - Gra ми Zycie 
https://en.wikipedia.org/wiki/John Horton Conwa 


ο Stephen Wolfram - New Kind of Science - 2002 
https://en.wikipedia.org/wiki/A New Kind of Science 
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Algorytm Ping-Pong 


e Inicjalizuj uktad 
ο Nałóż warunek początkowy (stan t=0) 


ο Dla każdej komórki ze stanu A (synchronicznie): 
° Oblicz regułę z sąsiedztwa 
e Zmień stan komórki 
* Zapisz jej wartość do stanu B 


° Skopiuj stan В do stanu А (lub zamień) 
ο Wizualizacja, obliczenia 
ο Powtórz pętlę 


Paski na zebrze? 
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Plamy - najprostszy model 


° Losowy stan poczatkowy р-0.5 


Algorytm 


ο Policz sume sasiadów komórki 
* Ustal jej nowa wartosc: 
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Gra w zycie 


Game of Life 


= | | * John Conway (1937-2020) 
io ΝΕ \ / | • Automat komórkowy 

s N —^ * Zywe/martwe komórki Тури 
© — | boolowskiego (0 lub 1) 
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-e- Publications (total) 
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INICJATYWA DOSKONAŁOŚCI 


Gra w życie 


John Conway (1937-2020) 


Więcej, niż gra w zycie 


"Larger than life" 


ο Rozszerzone sąsiedztwo 
ο Bardziej skomplikowane reguły przeżywania i narodzin 


ο Powstają większe struktury 


"Larger than life" 


ο Sasiedztwa: 


+ O 
EN Ep. 


e Cata sztuka polega na dobraniu odpowiedniej funkcji przezycia. 


Larger than life 


° R=5 ("bugs" - kwadratowe sasiedztwo) 


° 5 >= 34 && $ <= 58 -> komórka przeżywa (1-21) 
ο 5>=34 && s<=45 -> komórka sie rodzi (0 lub 1 -> 1) 


e Ciekawe źródła: 


https://www.youtube.com/c/slackermanz 
https://tatasz.github.io/compound_ca/ 
https://youtu.be/RZyy5bwXrn8 


Larger than life 


Wrażliwość па sąsiedztwo і szczegóły 
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Sasiedztwo kotowe, r=5, warunki: Sasiedztwo kwadratowe, r=5, warunki: 
przeżywa dla 34 <= $ <= 58, rodzi sie 34«s«45 przezywa dla 34 <- 5 <= 58, rodzi sie 34<=s<=45 


Glidery w "Larger than life" - 
https://youtu.be/wRGPZI7NcKI 


Multiple Neighborhood Cellular Automata 


* ММСА - uogölnienie automatow komórkowych 
“ Duza dowolnosc му doborze sasiedztw 
е Możliwe stany binarne lub ciągłe (0/1 lub 0-1) 


Zachowania odkryte νν MNCA 


° Tworzenie pojedynczych "osobników" 

ο Sity oddziatywania (przyciaganie / odpychanie) 

° Tworzenie silnych wiązań między obiektami 

° Ruch kolektywny (flocking) 

- Replikacja form 

- Oddziaływania elastyczne 

° Tworzenie kryształów 

e W niektórych wariantach spełnione zasady zachowania (пр. masy) 


(za https://slackermanz.com/understanding-multiple-neighborhood-cellular-automata/) 


Przyktad 1 MCNA, "cells" 


ο Dwa sąsiedztwa 
ο Sąsiedztwo 1: koło o promieniu 3 
e Sąsiedztwo 2: koło o promieniu 7 bez sąsiedztwa 1 


"Cells" 


if (avg >= 0.120 && avg <= 0.150) {res = 0; } 

if (avg >= 0.210 && avg <= 0.220) {res =1;} L Są εἰεοία Ἕως A 
if (avg >= 0.350 && avg <= 0.500) {res = 0; } 

if (avg >= 0.750 8.8. avg <= 0.850) {res = 0; } 


if (avg2 >= 0.100 8.8. avg2 <= 0.280) £ res=0;) ς sę امع‎ Xo 2, 
if (avg2 >= 0.430 && avg2 <= 0.550) { res = 1; } 


Algorytm 


1. Inicjalizuj tablice 
222. Sumuj wartości z sąsiedztwa 
3. Oblicz średnie avg/avg2 
4. Zastosuj reguły (umieranie, przeżywanie) 


5. Rysuj i powtarzaj od 2. 
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"Cells ММСА" 
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Przeglad: 


https: 


slackermanz.com/understanding-multiple-neighborhood-cellular-automata 
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Generalization of Conway 8 "Game of Life” to а 


continuous domain - SmoothLite 


Stephan Rafler 
Nürnberg, Germany 
frIndmr@web.de 


December δ. 2011 


Abstract 


We present what we argue is the generic generalization of Conway's 
"Game of Life" to a continuous domain. We describe the theoretical model 
and the explicit implementation on a computer. 


https://arxiv.org/abs/1111.1567 
- przeglad pracy 


SmoothLife 


* Model podobny do "Larger than Life" MNCA 
e Wartości 0-1 w komórkach (ciągłe) 
ο Całkowanie po dwóch przyległych otoczeniach min 


© 


Game of Life vs SmoothLife 


ο Wartość komórki - zywa/martwa 


ο Sasiedztwo 
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SmoothLife 


we can determine the filling of the cell or "inner filling" m by the integral 


1 
m = ዝተ И du f (X + t,t) 


\äl<r; 


and the neighborhood or "outer filling" n by the integral 
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where N and M are normalization factors such that the filling is between 


SmoothLife 


ο Funkcja przejscia ("transition function") postaci o(n,m) 


smoothness. For example we could define 


1 
1 + exp(-(x — а)4/а) 


01(2,0) = 


б2(2,04,0)--01(2,а)(1-01(2,9)) 


Ст(2,у,т) = 2(1 — o1(m,0.5)) + усі (т, 0.5) 


then we can define the transition function as 


s(n,m) = σο(η, ст (bi, 41, т), ст (bo, da, m)) (6) 


where birth and death intervals are given by [61, 65] and [41, d>] respectively. 


(uwaga: w praktyce lepiej zamienic sigma2 2 sigmam - o tym dalej) 


o1(x, a, а) 


1 


sigma1(x,0.5,0.1) σι(1,α) = 


sigma1( x, a, alpha) = 1.0 / (1.0 + exp(-(x - а) * 4 / alpha)) 
plot [0:1] 01(х,0.5,0.1) 


sigma2(x,0.5,0.6,0.1) —— σο(τ.α. Б) = σι(τ,α) (1 — oi(z,b)) 


sigma2(x,a,b,alpha)=sigma1(x,a,alpha)*(1-sigma1(x,b,alpha)) 
plot [0:1] 01(х,0.5, 0.6, 0.1) 


mit, у, т) = x(1— 01(m,0.5)) + ус (т, 0.5) 


ст(р, а, т, а) 


sigmam(0.257,0.365,x,0.01) 


sigmam(b,d,m,alpha) =b*(1-sigma1(m,0.5,alpha)) + d*sigma1(m,0.5,alpha) 
plot [0:1] [0:1] sigmam(0.257,0.365,x,0.01) 


s(n, m) = G2(N, Gm (01, di, m), Gm (22, do, m)) 


s(n,m) 


s(x,y,0.028,0.097,0.257,0.365,0.336,0.549) 


s(n,m,alphan,aloham,b1,b2,d1,d2)=sigma2(n,sigmam(b1,d1,m,alpham), 
sigmam(b2,d2,m,alpham),alphan) 


splot [0:1] [0:1] s(x,y,0.028,0.097,0.257,0.365,0.336,0.549) w pm3d 


s2(n,m) - inna postać, inne zachowanie modelu 


s2(x,y,0.028,0.097,0.257,0.365,0.336,0.549) 


5(х,у,0.028,0.097,0.257,0.365,0.336,0.549) ——— 


o2=om( o2(n, 01, b2, ап), с2(п, d1, d2, ап), m, am) 


s2(n,m,alphan,alpham,b1,b2,d1,d2)= 
sigmam(sigma2(n,b1,b2,alphan),sigma2(n,d1,d2,alphan), т, alpham) 


Nowy stan w SmoothLife 


new = old + dt * (2 s(n,m) - 1) 


Petny algorytm SmoothLife 


e |nicjalizuj mapę "życia" f(x,y) - wartości 0-1 
«Oblicz stopień zajętości sąsiedztwa 1 (m) 
«Oblicz stopień zajętości sąsiedztwa 2 (n) 
e Oblicz funkcje sigma1/2/m 
«Wyznacz nową wartość: 

* new = old + dt * (2 s(n,m) - 1) 
* Normalizuj "new" do wartości 0-1 


SmoothLife 


https://youtu.be/PUsWxzfuUmA 


Basic SmoothLife* 


https://www.shadertoy.com/view/sdKSzV 


52(х,у,0.028,0.097,0.257,0.365,0.336,0.549) 
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SmoothLifeS* - bariery aktywacji 


М<0.5 (martwa) N 
O st 
(М 20.25 && N<0.33) i о. 2» 
new = old + dt; ese 
else ~ 0.2 
new = old - dt; q 
N 0.57 O 


M>0.5 (Zywa) 
if(N>0.38 && М«0.56) 
new = old + dt; 
else 
new = old - dt; 


https://www.shadertoy.com/view/sdKSzV 


Simple SmoothLife 


e Komórki częściowo zywe/martwe 


ο Dwa sąsiedztwa 


ο Niecatkowite bariery dla przeżycia: 


Dodanie życia: (m«0.5, n>0.25, n<0.33) lub (m>0.5, n>0.35, п<0.51) 
Odjęcie życia: (m«0.5, n>0.25, n<0.33) lub (т<0.5, n>0.25, п<0.33) 


SmoothLifeS - algorytm 


1) Przygotowanie warunku poczatkowego (wartosci od 0 do 1) 
2) 
3) 

| 


4) Dodanie / odjecie wartosci zgodnie 7 barierami aktywacji 


Obliczenia gęstości żywych komórek w otoczeniu punktu m 


j.w. w otoczeniu n (torus 2D) 
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(kolorowanie wg M*N) 
https://youtu.be/6Wh5E6VerGA 


Stowo na koniec 


ο Dlaczego R1/R2=3? Np. R1=4, R2=12? 


- 08:91 9/9/5፻0 
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Wynik dla kwadratowego? (wstawic) 


N1/N2 Dla gry ми zycie 
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Analiza modelu SmoothLife 


ο Rozkład komórek w przestrzeni 

ο Klastrowanie faz 

* Obiekty zanurzone w środowisku (otoczenie, które je chroni) 
Jaka jest dynamika tych zmian? 
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"Diagram 
fazowy" 
N vs M 


"Diagram 
fazowy" 
N vs M 
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Koniec 


